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ABSTRACT 

Aims. We present a new upgraded version of ARES. The new version includes a series of interesting new features such as automatic 
radial velocity correction, a fully automatic continuum determination, and an estimation of the errors for the equivalent widths. 
Methods. The automatic correction of the radial velocity is achieved with a simple cross-correlation function, and the automatic 
continuum determination, as well as the estimation of the errors, relies on a new approach to evaluating the spectral noise at the 
continuum level. 

Results. ARES v2 is totally compatible with its predecessor. We show that the fully automatic continuum determination is con¬ 
sistent with the previous methods applied for this task. It also presents a significant improvement on its performance thanks to the 
implementation of a parallel computation using the OpenMP library. 
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1. Introduction 


One of the most commonly used spectral analysis methods is 
the line-by-line analysis, which in most of the cases requires a 
precise measurement of the strength of the spectral lines. This 
is normally done by measuring the equivalent-width (EW). EW 
methods have been used for decades in spectroscopy. 

In the past decade several automatic codes have emerged 
that aimed at releasing the astronomer from this tedious task 
(when using interactive routines), and, more importantly, to 
allow the precise and systematic measurement of equivalent 
widths. ARES (ISousa et alJ'2007 ) was one of the first codes to 
appear freely available to the community. ARES measures the 
equivalent widths of lines in an automatic way, trying to follow 
as closely as possible the same steps as those that would be taken 
by an astronomer when using interactive routines. Several other 
codes exist and appeared more recently to perform th e sam e 
task u sing different approaches: DAOSPEC dStetson & Pancinol 
DOOp, which is an automat ic wrapper for D AOSPEC 
(ICantat-Gaudin et aPl2014l) : TAME dKang & Leell2012l) . which 
follows the basic ideas from ARES very closely, is writ¬ 
ten in IDL, and presents some interactive facilities for the 
spectral analysis. Rencently, a more complete spectral anal¬ 
ysis code way presented that also computes EWs (iSPEC - 
Blanco-Cuaresma et alJl201~4h , 

Because it is so simple to use, ARES has been one 
of the most wi dely adopted codes for automatic EW mea¬ 
surements (e. g., Hekker_&Melendez l2007t iNeves et~aD [2009; 


task u 

2008 b 


ISozzetti et alJ 120091: iGhezzi et aL 20101: iRuchti et a~ 


2011 


* Automatic Routine for line Equivalent widths in stellar Spectra - 
ARES Webpage: http://www.astro.up.pt/~sousasag/ares/ 

** Based on observations made with ESO Telescopes at the La Silla 
Paranal Observatory under programme ID 075.D-0800(A) 


Sousa et al.1201 li Adibekvan et alJ2012a 

Tabernero et al. 

2012; 

Mortieretal. 2013^ Santos et al. 2013 

Monroe et al. 

2013; 

Maldonado et al. 2013; Jofre et al. 2014) 

ARES has also been 


used for astero seismology studies that require stellar charac¬ 
terization (e.g . Creevev et al.l 120121 : ISilva Aguirre et all 120131: 
Metcalfe et al. 2010) and is used i n large spectroscopic surveys 
such as the GAIA-ESO Survey dGilmore et al.l [20121) to ana¬ 
lyze the high-resolution spectra l data i n several working group s 
of the consortium (e.g. , ISousal d2014li : ISmilianic et al.l ( 20141) : 


lLanza fame et al. (2015)). The algorithms in ARES are also be¬ 
ing implemented in recipes for the very first data analysis soft¬ 
ware that will be available for a VLT instrument, ESPRESSO 
(Meg evand et al.l [201 2l iDi Marcantonio et al.l 120 1 2i iPene et all 
120 U 


We present the new improved version of the ARES code with 
special focus on the description of the new features that were in¬ 
troduced. A more detailed description of new features is avail¬ 
able on the ARES webpage. 


2. New features 

Correction for radial velocity: In the previous version, the in¬ 
put spectra needed to be fully corrected for radial velocity. Oth¬ 
erwise, ARES would not be able to correctly identify the spectral 
lines that were to be measured. In the new version, the user has 
two options that can be chosen to correct for the radial velocity 
in ARES: The first is to provide the value for the radial veloc¬ 
ity of the star in km/s. The second is to provide a mask with a 
list of lines at rest frame that will then be used by ARES to es¬ 
timate the radial velocity. To provide this new input to ARES, 
the new parameter rvmask needs to be included at the very end 
of the input file mine.opt. An example and a full description of 
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the method used for this new feature is provided on the ARES 
support web-pagcfl 

Automatic rejt estimation: The parameter rejt is the most 
important input parameter for the proper normalization and is 
crucial for the correct EW measurement, rejt strongly depends 
on the signal-to-noise ratio (S/N). Originally, this was a free pa¬ 
rameter set by the user because the visualization of the contin¬ 
uum position is highly subjective. This parameter continues to 
be free in ARES v2, and we still recommend this as the best op¬ 
tion for the user who desires to derive precise and accurate EWs. 
In some of our previous works we have presented som e sug ges- 
tions for the best value for this parameter. First, in iSousa et al.l 
(2008) we provided a discrete relation between the spectrum S/N 
and the recommended rejt value. This table gives only values for 
S/N > 100. ISousa et all ( 2011 1 later provided an empirical re¬ 
lation for the cases where the S/N is low (S/N < 100). These 
relations were derived for HARPS spectra, but according to our 
experience, they can be safely applied to other spectra whose 
S/N was computed in a similar way. If the S/N is not available, 
we recommend either adjusting this parameter individually for 
each spectra or finding a way to first estimate the S/N for the 
spectra and then derive the proper relation with the rejt parame¬ 
ter (for an example of such an exercise see lMortier et all2013l) . 

The new version of ARES has an additional option where 
the user can specify spectral regions in the spectra, which should 
ideally be clear of spectral lines. ARES will internally estimate 
the S/N and relate it with a proper rejt value. For a detailed de¬ 
scription of this new option, see Sect. [3] 

rejt dependence on wavelength: The S/N can change signif¬ 
icantly over the wavelength range for some spectra. With this in 
mind, we have implemented a way to change the rejt parameter 
for different wavelengths. To use this option, the rejt parame¬ 
ter has to be defined with the value -2. ARES will then read a 
file with the name “lambda_rejt.opt\ which should be an ASCII 
file with two columns (wavelength and rejt separated by spaces). 
ARES will read the values in the file and then interpolate the rejt 
value for each wavelength and will adjust the continuum accord¬ 
ingly for each spectral region. 

EW error estimation: One of the caveats in the original 
code was that it lacked an error estimation or the flagging to 
identify problematic measurements. This new feature will cer¬ 
tainly be helpful, especially in spectral analyses where only very 
few lines are available. We have implemented an efficient esti¬ 
mation of the error for the EWs. For more details on this new 
task see Sect. [4] 

Read ASCII spectral data: Until now, ARES was only able 
to read specific FITS spectral data. Since there are several spec¬ 
tra available in public archives that follow different formats (e.g., 
the ESPADONS or NARVAL archives), the new ARES code can 
now read ASCII data. If a FITS file is not provided, ARES will 
try to read the file as ASCII data. For this to work, the data should 
be available in two distinct columns (wavelength and flux sepa¬ 
rated by spaces). 

Cleaning any zero spectral flux: One of the first tasks in 
ARES is to determine whether the input spectra have zero values 
in the flux. This usually appears in spectral data from spectro¬ 
graphs that have a gap in the spectral coverage. In many cases, 
the standard procedure for these cases is to assign zero values 
for these spectral regions. Unfortunately, zero values in the flux 


1 http://www.astro.up.pt/ sousasag/ares/support.html 


can be problematic for the spectral normalization, and these re¬ 
gions should be avoided to prevent the code from crashing. To 
overcome this problem, ARES scans the spectral data and re¬ 
places them for any zero value by unity, which avoids unex¬ 
pected crashes. 

Parallelization of the code with OpenMp: ARES is now 
implemented to use OpenMp. This is a code parallelization li- 
brar\Q that simplifies the interface with multicore machines. On 
the ARES webpage we provide some benchmark examples that 
were used to compare and test the performance of the new ver¬ 
sion with the older one. To summarize, these tests improve the 
performance from ~ 60% for K stars up to ~ 250% for the F and 
G stars. 

All these features are now implemented in ARES in a way 
that the code remains fully compatible with the previous version. 
The user is free to choose or discard any of the new features ex¬ 
cept for the errors, which will always appear in the new ARES 
output without any significant loss of performance; quite to the 
contrary, performance will improve if the code is used in a mul¬ 
ticore machine. 

3. Automatic rejt 

The identification of the continuum position is a very difficult 
and subjective task in these spectral analyses. The rejt parameter 
plays a fundamental role in ARES for this purpose (see Fig. 2 
and Sect. 3.1 of lSousa et al.l (20071) for a complete description of 
this parameter). Because of that, this parameter still remains ad¬ 
justable for the user. One of the advantages of automatic codes 
such as ARES is that when the parameter for the continuum is 
defined, the code will measure all lines in the spectra in a sys¬ 
tematic way, in this way avoiding errors that are caused by the 
subjectivity that is common in the use of interactive routines. 

To keep the EW measurements between different spectra ho¬ 
mogeneous, the best option is to use an empirical relation be¬ 
tween the spectral S/N and the rejt parameter. This ensures the 
homogenization of the measurements for the analysis of sev¬ 
eral similar spectra that were observed in the same spectroscopic 
configurations. 

This solution was unsatisfactory, therefore we added the pos¬ 
sibility of computing an automatic rejt parameter. This is a diffi¬ 
cult task, and some caution should be taken. This option can be 
accessed by defining the rejt parameter in a different way. For ex¬ 
ample, setting r<yY=3;5764,5766,6047,6052,6068,6076 will tell 
ARES to evaluate the noise in three different spectral ranges, 
the first number (3) defines the number of ranges, and each next 
pair of numbers represents each wavelength region. These spec¬ 
tral regions should be selected carefully and avoid strong lines. 
The ideal case will be without any lines at all. The three spectral 
regions given in this example are regions that have no lines, or 
at most very week lines, for solar-type stars, and therefore we 
recommend them. 

3.1. Estimating the signal-to-noise ratio 

To obtain an automatically determined rejt parameter, we start 
with estimating the S/N. Figure Q] is used here to explain the 
steps of estimating the S/N: 

- The first step is extracting the spectral region. The original 

spectrum is represented in blue for all panels, and the region 

is marked by the gray area in the panels. 

2 http://www.openmp.org 
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Fig. 1. Evaluation of the noise in the three spectral regions recommended to automatically estimate the rejt parameter (see text for the details). 


- The second step is to smooth the original spectra. This is 
achieved through a boxing car average. This allows filtering 
the noise at higher frequencies (in wavelength) (for more de¬ 
tails see Sect. 3.3 of ISousa et al.l[20071) . The value for the 
smoothing in this procedure was hard-coded in ARES with a 
value of eight pixels. The smooth spectra are represented in 
green shifted down by a constant for clarity. The idea is that 
these spectra do not represent any lower frequency signal, 
where the spectral lines are included. 

- Next we subtract from the original spectrum this smoothed 
spectrum that is free from high-frequency noise and place it 
at the spectrum level (using the average of the flux for each 
spectral region). This is represented with the red lines, again 
slightly shifted by a constant for clarity. These data in red 
represent the high-frequency noise that was used to estimate 
the S/N for each spectral region. 

- The S/N computed for each spectral region is simply its mean 
flux < F > divided by the standard deviation of the sub¬ 
tracted spectra cr F (spectra in red): S/N — 

- The final adopted S/N value is derived by taking the median 
of the three values derived for this example. The choice of 
the median is more preferable than the average since it is not 
affected by the outlying points (e.g., a cosmic-ray hit). 

Figure [I] shows the estimation of the S/N for the three spec¬ 
tral ranges given as example. Here we used a HARPS spectrum 
from the archive for the star Alpha Cen A. This spectrum was 
observed on 2005-04-23, 09:57:42 UT. A S/N of ~ 275 is indi¬ 
cated in the header of this spectrum. This value is compatible 
with the median of the values presented in the figure (256.29). 

3.2. Auto-rejt and ARES local continuum determination 

The simple relation between rejt, the S/N, and the standard de¬ 
viation of the noise is re jt - 1 - = 1 - cr F , using < F >- 1 

for the normalized flux at the continuum. 

The auto -rejt option should be used with caution since the 
spectral ranges selected may be contaminated with unexpected 
spectral features, such as cosmic rays or telluric lines. The S/N 
can also be directly included in the rejt parameter. In this case, 
the S/N is not estimated and directly uses this relation. This is an 
excellent option if the user has a prior knowledge on the S/N of 
the spectra. 

To test the new feature, we present in Fig. [2] a comparison 
of EWs determined in two different ways. Here we used seven 



Fig. 2. Comparison of the equivalent widths measured using the recom¬ 
mendation from lSousa et all d2008t) for the rejt parameter and using an 
automatic value. 

spectra with different S/N, ranging from -70 to -2000, measur¬ 
ing ~ 380 lines for each of them. The goal is to show the dif¬ 
ferences betwen using the new auto -reU option with the option 
of using the recommended values from lSousa et al.l (l2008t) . This 
comparison shows that there is a small mean offset of - 0.7 mA 
with a still significant dispersion of - 4.4 mAMost of the large 
spread is due to the spectra with lower S/N. The panel inside the 
figure shows a zoomed part of the comparison where it is clear 
that the high percentage (~ 96%) of the lines are indeed close to 
the identity line. 

We would like to recall that the importance of these auto¬ 
matic codes is their systematic performance. Regardless of the 
option selected in ARES or the value adopted for the rejt pa¬ 
rameter, the code will always behave systematically. This is very 
important for instance for a differential analysis, for which it is 
crucial that exactly the same method is used in the calibration 
and the spectral analysis. Another good example where the sys¬ 
tematic measurements are very important is the use of line ratios 
to der ive spectral indexes or even parameters such as the temper¬ 
ature (iKovtvukh et al.l[2003USousa et alJl2010Ll2012h . 

4. EW error estimation 

There are many sources contributing to the EW error. These in¬ 
clude the determination of the continuum position, the S/N of 
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the data, blended lines, or any unexpected spectral feature, such 
as cosmic-ray hits. 

An easy and fast solution was implemented in ARES to es¬ 
timate the error on the equivalent width for each line. The error 
estimation is directly based on the fit that ARES performs on the 
spectral lines. Recalling this process, ARES fits the spectral lines 
with Gaussian profiles, using a non-linear square fit algorithm 
(gsl_multifit_fdfsolver_lmder - defined in the Gnu Scientific li¬ 
brary, GSL, as a robust and efficient version of the Levenberg- 
Marquardt algorithrrQ). To use this algorithm, the Gaussian pro¬ 
file and the respective parameter derivatives (the so-called Jaco¬ 
bian matrix) was implemented in the code. The Gaussian profile 
implemented in ARES is written as 


G(t) = sM = Yj ^' ex P(- A <( f - G)). 


( 1 ) 


where di is the depth of the ith line, c, its center, and A, is related 
with the Gaussian width such as A,- = 1 /(2cr 2 ), where cr is the 
Gaussian width (also known as the standard deviation). This spe¬ 
cific representation of the Gaussian profile was chosen to closely 
follow the GSL example for the fitting of a single Gaussian with 
the non-linear square fit algorithm. The difference to the previ¬ 
ous ARES version is that the fit now considers an error for the 
spectral data. This is then used to estimate the errors on each 
parameter for the Gaussian fit. Note that although some reduced 
spectra in archives do have the error on the flux for each pixel, 
many do not. Because of this, and also to ensure consistency 
with the previous version of ARES, we still only require the in¬ 
put of the wavelength and the respective flux for each pixel. The 
error adopted for each pixel is taken as the same (every pixel 
will have the same weight), and its value strongly depends on 
the S/N, which is derived using exactly the same procedure as 
described in the previous section ( cr F ). If the rejt parameter is 
defined by the user, <x F will be derived from its relation with the 
S/N. 

With the uncertainty given for each of the fitted parameters, 
we can use it to estimate the error of a measurement. First, the 
equivalent width in ARES is computed directly from the area of 
the fitted Gaussian profile: 


EW; = 


■-r 


gi(t)dt = diCTi '{in = di J -L 

V 


( 2 ) 


where i represents the Gaussian that fits the line to be measure. 
Now we can use the uncertainties derived from the fit on the 
line depth (Adi) and on the line width (AA,), and from the error 
propagation equation we derive 


A EWi = EWi 


Adi 

di 


1 AA/ 

2 A“ 


(3) 


This error estimation includes the error on the continuum posi¬ 
tion (through the use of cr F ) and also part of the error introduced 
by blended lines that are responsible for the larger errors of the 
Gaussian profile parameters. 

The main advantage of this new feature is that it is possible 
to immediately identify problematic EW measurements. This is 
very useful for spectral analysis, which relies on very few lines. 
For example, the analysis of elements such as Na, Mg, Al, Sc, 
and M n are typically based on ly on two to four lines (per ele¬ 
ment) (Adibek van et aljfeoi 2bl) . The studies based on EW mea¬ 
surements for these elements using a limited line list will cer¬ 
tainly benefit from evaluating the uncertaities for each line. 


5. Summary 

We presented a new version of the ARES code. This version 
includes new features. The most important are the automatic 
correction for radial velocity, the new option for automatically 
defining the continuum position, and the error estimation for the 
equivalent width measurements. All these features will certainly 
be useful for the users of ARES, especially for the analyses per¬ 
formed on large sets of spectral data. 

The errors derived in ARES are good enough for a prelimi¬ 
nary detection of poor measurements. This is a new capability, 
which can be used to pre-select the best-fit lines. We still recom¬ 
mend the user to determine whether there are any outliers before 
using them for further analysis. 

A more detailed description of the new features can be found 
on the ARES webpage. A proper understanding of these new 
features is fundamental for a proper judgment of the best option 
for each spectral analysis project. 

Finally, we would like to emphasize again that ARES v2 re¬ 
mains fully compatible with the previous version and will also 
work faster, thanks to the use of parallel processing for multicore 
processors. 
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